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We study the structure of eigenstates in two-body interac- 
tion random matrix ensembles and find significant deviations 
from random matrix theory expectations. The deviations are 
most prominent in the tails of the spectral density and in- 
dicate localization of the eigenstates in Fock space. Using 
ideas related to scar theory we derive an analytical formula 
that relates fluctuations in wave function intensities to fluctu- 
ations of the two-body interaction matrix elements. Numeri- 
cal results for many-body fermion systems agree well with the 
theoretical predictions. 
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Random matrix theory (RMT) has become a powerful 
tool for describing statistical properties of wave functions 
and energy levels in complex quantum systems j^j^ . The 
use of a two-body random matrix ensemble (TBRE) [^|^] 
is of particular interest for many-body systems since clas- 
sical RMT implies the presence of fc-body forces {k > 2) 
and gives the unphysical semicircle as the spectral den- 
sity. The TBRE displays the same spectral fluctuations 
as classical RMT while its Gaussian spectral density 
agrees well with nuclear shell model calculations. The 
situation is not so clear for the structure of wave func- 
tions in TBRE. Recent results show that ground states of 
shell model Hamiltonians with two-body random inter- 
actions favor certain quantum numbers and thereby 
differ considerably from RMT expectations. 

In this letter we examine the structure of wave func- 
tions in TBRE, and compare numerical results with the- 
oretical predictions. Such a study is not only interesting 
on its own but is also motivated by the ongoing impor- 
tance of TBRE for nuclear |^,|| and mesoscopic physics 
. We recall that deviations from RMT indicate some 
degree of wave function non-ergodicity and are related 
to phenomena like Fock space localization in many-body 
systems P,|l0|, and scars of periodic orbits |ll| , [l^ or in- 
variant manifolds in classically chaotic systems. 

To quantify the degree of localization of a given wave 
function, it is useful to introduce the notion of an in- 
verse participation ratio (IPR). Thus, let D be the to- 
tal dimension of the relevant Fock subspace, let \b) label 
the single-particle basis states (Fock states), and let \a) 
represent the eigenstates of the Hamiltonian. Then the 
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are the squares of the expansion coefRcients, and have 
mean value 1/D. The IPR of eigenstate |a) is defined as 
the first nontrivial moment of the intensity distribution, 
namely the ratio of the mean squared Pat to the square 
of the mean: 
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The IPR measures the inverse fraction of Fock states that 
participate in building up the full wave function \a), i.e. 
IPRa = 1 for a wave function that has equal overlaps 
Pab with all basis states, and IPR^ = D for the other 
extreme of a wave function that is composed entirely 
of one basis state. While complete information about 
wave function ergodicity is contained in the full distri- 
bution of intensities Pabi the IPR serves as a very useful 
one-number measure of the degree of Fock-space localiza- 
tion. In RMT, the Pab are given (in the large-D limit) by 
squares of Gaussian random variables, in accordance with 
the Porter-Thomas distribution, leading to IPRrmt — 3 
for real wave functions. For finite D, the IPR is slightly 
below its asymptotic value (with the deviation falling off 
as 1/-D), but is still uniform over the entire spectrum. 

A very simple two-body interaction model, however, 
already displays behavior which is qualitatively different 
from this naive RMT expectation. Let the Hamiltonian 
be given by 
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where the single-particle indices i, j, k, I run from 1 to 
M (the number of available single-particle states) and 
the VijM are Gaussian random variables with unit vari- 
ance. The operators and aj create and annihilate a 
fermion in the single-particle state labeled by j, respec- 
tively, and obey the usual anti-commutation rules. The 
dimension of the iV-particle Fock subspace is given by 
D — (^) . We notice that this model contains no explicit 
one-body terms and therefore does not display Anderson- 
type localization effects. On physical grounds, we assume 
that the two-body interaction is generated by a real sym- 
metric potential, U(ri,r2) = U(r2,ri) and choose real 
single-particle wave functions. These conditions lead to 
the constraints Vikji — Vijki — Vjuk = Vijki and reduce 



1 



the number of independent variables describing a given 
realization (Eq. |^) of the ensemble. While such correla- 
tions between matrix elements will determine factors of 2 
in the calculation below, they do not qualitatively affect 
our results. The choice of bosons rather than fermions 
also does not qualitatively change the localization behav- 
ior. 

Fig. |l](top) shows a smoothed ensemble-averaged spec- 
tral density 



p{E) ^ Tr S{E - H) = J2SiE ~ E^) , 
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for M — li and = 6, while in Fig. |l|(bottom) we plot 
IPR/IPRrmt — 1 as a function of energy. Obviously, 
we observe strong deviations from RMT behavior at the 
edges of the spectrum. This implies that wave function 
intensities do not have a Porter-Thomas distribution. 




FIG. 1. Spectrum and IPR for A'^ = 6 fermions distributed 
over M = 13 orbitals. Data is averaged over an ensemble 
of 20 systems of the form (Eq. Top: Normalized spec- 
trum p{E)/D, Eq. ^ (solid curve), Gaussian prediction, Eq. |l^ 
(dashed curve), and modified Gaussian, Eq. hq (dotted curveyT 
Bottom: IPR/IPRrmt — 1: data (solid curve), simple theory, 
Eq. |l^ (dashed curve) and theory based on modified Gaussian 
form (dotted curve). 

To understand this surprising result, we may adapt 
a formalism previously used successfully to understand 
wave function scars and other types of anomalous quan- 
tum localization behavior in single-particle systems. Let 
Pb{E) be the local density of states (strength function) 
of the basis state b: 



Pb{E) = {h\5{E - H)\h) = J2 PcbSiE - E^) ■ (5) 



it is given by the Fourier transform of the autocorrelation 
function 



A,{t) = (6|e-''^*|6) 
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Let us assume that A{t) displays only two different 
time scales: the initial decay time Tdocay of the Fock 
state 1 6) due to interactions, and the Heisenberg time 
Tji ^ ^Tdocay (i-G. ?i ovcr the mean level spacing) at 
which individual eigenlevels are resolved. Following the 
initial decay, random long-time recurrences in A},(t) can 
be shown to be convolved with the short-time behav- 
ior |l|,|l2|. In the energy domain this produces random 
oscillations fb{E) multiplying a smooth envelope pf^{E) 
given by the short-time dynamics: 



Here 



PbiE)^ME)priE). 



D 



(7) 



(8) 



while fb{E) is a fluctuating function with mean value of 
unity: 
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The Tab are random variables with mean value one. 
Then, substituting into Eq. |^, we have the individual 
wave function intensities given by 



Eab — ^ai 
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Both Tab and have 6-dependent fluctuations, which 
are uncorrelated under our assumption of separated time 
scales. Then using Eq. |o| we can express the IPR (Eq. H) 

as {5pr{E)^pT[E)~{pr{E))) 



(Tab)' {priEa))' 

{SpriEc.?} 



= IPRr 



MT 



(priE^))' 



(11) 



where all averages are over the Fock basis index 6, i.e. 

In the limit of many particles (and many holes) N, M— 
N ^ 1 the spectrum approaches a Gaussian shape 16 
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where E^ 



D~^Tr is given by the mean sum of 
squares of matrix elements in any given row of the Hamil- 
tonian. The same arguments lead to the vanishing of 
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higher-order cumulants for the individual strength func- 
tions p^J^{E), so each of these should also have a Gaus- 
sian shape, but with centroid ct, = iJ^f, and variance 
Vb = J2b' ■ ^° these fluctuations in the centroids 
and widths as one goes through the different basis states 
b, we have 
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where in the second line we have used the Gaussian form 
ofEq. Substituting into Eq. |ll|, we obtain the general 
form: 
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valid of course not only for two-body interactions but 
for any ergodic Hamiltonian with a Gaussian density of 
states. The quantities ((Jcf,)^) and {{Svb)'^) do depend 
on the parameters of the model, but remarkably the IPR 
behavior is always given to leading order by a quartic 
polynomial in energy. From this point of view, it is easy 
to understand the enhancement in fluctuations near the 
edge of the spectrum (i.e. for large \E\): two Gaussians 
differing only slightly in their centroid or width may look 
almost identical in the bulk, but the relative difference 
increases dramatically as one moves into the tail of the 
spectrum. 

For our simple model (Eq. |^) one may easily compute 
the coefficients in Eq. |l^. The squared width (vb) — Eq of 
the full spectrum is given by the sum of squares of entries 
in one row of the Hamiltonian and equals the number of 
independent terms in Eq. |^ that couple any given Fock 
state to other Fock states, times the mean squared value 
(= 1) of each such term. A simple counting argument 

then shows 
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cause all the contributions are independent variables 
of mean 1 and variance 2, the variance in the sum is given 
by twice the number of contributions: ((^Wb)^) — lE^. 
Finally, the variance in the centroid is given by the num- 
ber of terms in the Hamiltonian (Eq. ^ which contribute 
to each diagonal element in the Fock basis, namely 
{[bcbf) = (^). Then 
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where in the last line we have taken the dilute many- 
particle limit M ^ N 1. The fluctuations are al- 
ways strongest near the edge of the spectrum; in par- 
ticular close to the ground state (E'g^ « 2El In D « 
2ElN\n{M/N)), we find 
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Let us compare the IPR prediction of Eq. with the 
numerical data presented in Fig. |l|(bottom). While the 
prediction qualitatively describes the correct trend of the 
IPR as a function of energy, it is not in quantitative agree- 
ment with the data and fails by as much as a factor of 
2 at the edge of the spectrum. This is not very sur- 
prising, given that the actual behavior of the spectral 
density p^'^{E) is far from Gaussian for the parameters 
we have chosen. Although it is known that the spec- 
trum approaches a Gaussian form in the many-particle 
limit | [T5| , |l6t , the number of particles that can be simu- 
lated numerically is not nearly large enough for the Gaus- 
sian form to be a good quantitative approximation, par- 
ticularly in the tail (see Fig. 0(top)). A good way to 
measure deviations from the Gaussian form is to com- 
pute the fourth cumulant divided by the square of the 
second cumulant E^: for a Gaussian this quantity van- 
ishes while for our system it equals —0.7, actually closer 
to the semicircle value of —1. 

To obtain quantitatively valid predictions, we must 
correct for deviations from the Gaussian shape. The 
ansatz 
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allows for nonzero higher cumulants while keeping the 
number of states and the width E^ fixed. The parameter 
e is obtained by a least squares fit of the numerical data 
to this form. A derivation starting from the first line of 



Eq. 13 yields an IPR that depends on the energy through 



a rational function instead of the quartic polynomial of 
Eq. |l5|. As we can see in Fig. |l](bottom), this leads to 
a surprisingly good quantitative prediction for the IPR 
behavior, given that we are applying perturbation theory 
around a Gaussian shape for a spectrum that in reality 
is very far from the Gaussian limit. 

Gaussian spectral densities have been observed in nu- 
clear shell model calculations with realistic [0 or ran- 
dom two-body interactions This is mainly due to 
the presence of spin and isospin degrees of freedom, which 
cause important correlations between Hamiltonian ma- 
trix elements. Let us therefore consider adding spin to 
the Hamiltonian (Eq. ||), so that a] and aj now create and 
annihilate a fermion in the single-particle state labeled by 
j = (rij, Sj), with rij and Sj — ±1/2 denoting the orbital 
and spin quantum numbers. In what follows we consider 
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N fermions in a shell of fixed total spin S = — ^■ 

We assume that the random matrix elements Vijki de- 
pend only on the orbital quantum numbers; this reduces 
the number of independent matrix elements considerably. 
The (sparse) two-body interaction matrix is constructed 
using a code similar to the one described in Rcf. . 

The spectral density of this TBRE agrees very well 
with a Gaussian, and the IPR is predicted by Eq. 
However, the determination of the quantities {{Scb)"^) and 
(((5wf,)^) is more difficult since the spin degree of freedom 
makes the required counting of matrix elements a non- 
trivial task. Alternatively, one may obtain these quanti- 
ties directly from the numerically generated Hamiltoni- 
ans. Fig. || shows that the numerical data and theoretical 
results are in good agreement. This confirms the validity 
of the simple theory derived in this letter. 



two-body random ensemble exhibit fluctuations that de- 
viate strongly from random matrix theory predictions, 
while good agreement is obtained in the bulk of the spec- 
trum. The numerical results agree well with the theoreti- 
cal prediction that is derived from arguments used in scar 
theory. In particular, we have presented a simple formula 
that relates fluctuations of the wave function intensities 
to fluctuations of the two-body matrix elements of the 
Hamiltonian. 

We thank George F. Bertsch for suggesting this study 
and for several very useful discussions. This research 
was supported by the DOE under Grant DE-FG-06- 
90ER40561. 




FIG. 2. Spectrum and IPR for A'' = 6 fermions with spin 
distributed over M = 7 orbitals. Data is averaged over an 
ensemble of 20 systems of the form (Eq. |^. Top: Normalized 
spectrum p(E)/D, Eq. ^ (solid curve) and Gaussian predic- 
tion, Eq. |l2| (dashed curve) . Bottom: IPR/IPRrmt — 1: data 
(solid curve) and theory, Eq. ^ (dashed curve). 

We note that our results are in qualitative agreement 
with nuclear shell model calculations using realistic inter- 
actions 19 -2^ and with calculations in atomic physics 
p2| . Our calculation shows that deviations from ergod- 
icity near the ground state do not require the presence of 
one-body terms in the Hamiltonian, though such terms 
may of course enhance the degree of localization. 

In summary we have studied numerically and analyti- 
cally the wave function structure in many-body fermion 
systems with random two-body interactions. Near the 
edge of the spectrum, wave function intensities of this 
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